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Abstract. We explore the possibility that, in the vicinity of an OB star, a prestellar core which would otherwise 
have formed an intermediate or low-mass star may form a free-floating brown dwarf or planetary-mass object, be- 
cause the outer layers of the core are eroded by the ionizing radiation from the OB star before they can accrete 
onto the protostar at the centre of the core. The masses of objects formed in this way are given approximately by ~ 
0.010M G (flj/0.3 kms~') 6 (A/L yc /10 50 s~')~ 1/3 (n /l0 3 crrr 3 )~ 1/3 , where ci\ is the isothermal sound speed in the neutral gas of 
the core, A/Ly C is the rate of emission of Lyman continuum photons from the OB star (or stars), and no is the number-density 
of protons in the HII region surrounding the core. We conclude that the formation of low-mass objects by this mechanism 
should be quite routine, because the mechanism operates over a wide range of conditions (10 50 s~' S yV Lyc S 10 52 s~', 
10cm -3 5: no 5= 10 5 cirT 3 , 0.2 km s" 1 52 ai 5= 0.6 km s" 1 ) and is very effective. However, it is also a rather wasteful way of 
forming low-mass objects, in the sense that it requires a relatively massive initial core to form a single low-mass object. The 
effectiveness of photo-erosion also implies that that any intermediate-mass protostars which have formed in the vicinity of a 
group of OB stars must already have been well on the way to formation before the OB stars switched on their ionizing radiation; 
otherwise these protostars would have been stripped down to extremely low mass. 
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. | 1- Introduction (2004a,b). Here we explore an alternative possibility, namely 

■ that in the central regions of clusters like the Orion Nebula 

s_j ■ The existence of brown dwarves was first proposed on theo- ™ , ■ ,. ■ f f 

H F F Cluster, isolated brown dwarves can form from pre-existing 

03 retical grounds by Kumar (1963a,b) and Hayashi & Nakano . , • , ■ , u u . . • • ■ 

& 3 ' 1 prestellar cores which are eroded by the ambient ionizing ra- 

(1963), but more than three decades passed before observa- ,. t . c , , c . ., , c c ,. c c 

' r diation held. Similar explanations for the formation of free- 

tions were able to establish their existence unambiguously a .. , . , , jut. 

b 1 floating low-mass objects have been proposed by Zapatero 

(Rebolo, Zapatero Osorio & Martin, 1995; Nakaiima et al., ^ . , , /OAAAN . v , . , 

x ' ±- Osorio et al. (2000) and Kroupa (2002), but they do not appear 

1995; Oppenheimer et al., 1995). Now they are routinely ob- tU u i a ■ a * u * 

FF - J to have been explored in detail before now. 
served (e.g. McCaughrean et al., 1995; Luhman et al., 1998; 



Wilking, Greene & Meyer, 1999; Luhman & Rieke, 1999 
Lucas & Roche 2000; Martin et al. 2000; Luhman et al„ 2000 
Bejar et al. 2001; Martin et al., 2001; Wilking et al, 2002 



If a dense neutral core is immersed in an HII region, the 
ambient Lyman continuum radiation drives an ionization front 



McCaughrean et al., 2002). It appears that a significant fraction lnto ±s core - The ionization front is normally preceded by a 
of brown dwarves in clusters are not companions to more mas- compression wave, and this may render the core gravitation- 
sive stars, but rather they are free-floating. Reipurth & Clarke all y unstable. At the same time, the flow of ionized gas off 
(2001) have suggested that this is because brown dwarves are the core reduces lts mass and ln J ects matenal and klnetlc en " 
ejected from nascent multiple systems before they can acquire ^ lnto the surrounding HII region. This situation was ana- 
sufficient mass to become Main Sequence stars, and this sug- ^ b ? D y son ( 1968 ; see also Kahn ' 1969 ) who ^S^d that 
gestion is supported by the numerical simulations of Bate et 11 could ex P lain the broadening and splitting of emission lines 

from the Orion Nebula. He also considered the circumstances 



al. (2002), Delgado Donate et al. (2003) and Goodwin et al. 



under which such a core might become gravitationally unsta- 
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Bertoldi (1989) and Bertoldi & McKee (1990) have de- 
veloped semi-analytic models for larger cometary globules 
created by the interaction of ionizing radiation with a pre- 
existing cloud on the edge of an HII region. They use these 
semi-analytic models to evaluate the acceleration of the glob- 
ule (due to the rocket effect; cf. Kahn, 1954; Oort & Spitzer 
1955) and its gravitational stability. Their results agree broadly 
with the numerical simulations of Sandford et al. (1982) and 
Lefioch & Lazareff (1994). More recently Kessel-Deynet & 
Burkert (2002) have modelled the radiation driven implosion 
of a pre-existing large cloud having internal structure. They 
use Smoothed Particle Hydrodynamics to handle the three- 
dimensional, self-gravitating gas dynamics, and find that small- 
scale internal structure delays, and may even inhibit, gravita- 
tional collapse. 

However, in this paper we are concerned - like Dyson and 
Kahn - with small-scale neutral cores immersed well within an 
HII region, in the vicinity of the exciting stars, rather than large 
cometary globules at the edge of an HII region. In particular, 
we wish to evaluate the competition between erosion of mate- 
rial from the outer layers of such a core (due to the ionization 
front which eats into the core), and the formation and growth of 
a central protostar (triggered by the compression wave driven 
into the core ahead of the ionization front). 

2. Overview of Model 

In order to keep free parameters to a minimum in this ex- 
ploratory analysis, we make a number of simplifying assump- 
tions. We assume that there exists a prestellar core in equilib- 
rium, and that a nearby massive star (or group thereof) then 
switches on its ionizing radiation instantaneously, and ionizes 
the surroundings of the prestellar core, on a time-scale much 
shorter than the sound-crossing time for the core (we justify 
this assumption later). The resulting increase in external pres- 
sure triggers the collapse of the core, from the outside in, and 
at the same time an ionization front starts to eat into the core, 
thereby eroding it from the outside, as it contracts. The final 
stellar mass is therefore determined by a competition between 
the rate of accretion onto the central protostar, and the rate of 
erosion at the boundary. Erosion ceases to be effective when 
the ionization front has eaten so far into the collapsing core 
that the newly ionized gas flowing off the ionization front is - 
despite its outward kinetic energy - still gravitational bound to 
the protostar. 

The evolution consists of three phases, starting from the 
instant t — when the surroundings of the prestellar core 
first become ionized. At this juncture a compression wave is 
driven into the core setting up a uniform, subsonic inward ve- 
locity field. When this wave impinges on the centre at time 
h, the first phase terminates and a central protostar is formed, 
which then grows steadily by accretion. At the same time an 
expansion wave is reflected and propagates outwards leaving 
an approximately freefall velocity field in its wake. At time t 2 
this outward propagating expansion wave encounters the in- 
ward propagating ionization front, and the second phase termi- 
nates. Thereafter the ionization front is eroding material which 
is falling ever faster inwards. Eventually, at time t 3 , the ioniza- 



tion front encounters material which is falling inwards so fast 
that it cannot be unbound by the ionization front, and the third 
phase terminates. The final stellar mass is the mass interior to 
the ionization front at time t 3 . 

In Section 3 we derive an expression for the diffuse ionizing 
radiation field in an HII region, and in Section 4 we formulate 
the speed with which the ionization front eats into the core. In 
Section 5, we describe the initial prestellar core and formulate 
the time t\ it takes for the compression wave to propagate to 
its centre. In Section 6 we describe the modified density and 
velocity field set up in the inner core by the outward propa- 
gating expansion wave, and formulate the time t 2 at which the 
outward propagating expansion wave meets the inward propa- 
gating ionization front. In Section 7 we formulate the time t 3 
at which the ionization front has eaten so far into the modified 
density field inside the expansion wave that it has encountered 
material which is infalling too fast to be unbound by the act of 
ionization. The material interior to the ionization front at time 
t 3 is presumed to constitute the final star. We can therefore ob- 
tain an expression for its mass. 



3. The diffuse ionizing flux 

We shall assume that the HII Region which engulfs the pre- 
existing prestellar core is excited by a compact cluster of OB 
stars which emits Lyman continuum (i.e. hydrogen ionizing) 
photons at a constant rate N Lyc . We shall also assume that most 
of the surrounding space is occupied by a uniform gas in which 
the density of hydrogen nuclei in all forms is n . The radius of 
the initial HII Region is therefore 



^HII — 



3M 



Lyc 



4 n a, nl 



1/3 



(1) 



where a„ 2 x 10~ 13 cm 3 s _1 is the recombination coefficient 
into excited states only; in other words, we are using the On- 
The-Spot Approximation. 

If we neglect the shadows cast by other prestellar cores and 
attenuation by dust inside the HII region, the radial - or direct - 
number- flux of Lyman continuum photons A^ lid (i.e. the num- 
ber of ionizing photons crossing unit area in unit time) varies 
with distance R from the OB cluster as 



"Lyc 



Lyc 



4nR 2 



(2) 



where the first term on the righthand side represents geometric 
dilution and the second term represents the absorptions which 
maintain ionization against recombination. Hence the mean 
value of N^ a , averaged over the volume of the HII region, 
is 



«c ial ) 



An 



(3) 



The volume number-emissivity for the diffuse component 
of the Lyman continuum radiation field (i.e. the number of 
ionizing photons emitted from unit volume in unit time into 
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unit solid angle, due to recombinations straight into the ground 
state) is 



'Lyc 



4n 



(4) 



where or, - 10 cm s is the recombination coefficient into 
the ground state only. The mean diffuse flux is therefore of or- 
der 



(n: 



diffuse \ 
Lyc / 



4a t 



f 3 A/ljc af nl » 



4n 



(5) 



Comparing Eqns. (0 and (|5}, we see that the mean diffuse flux 
is smaller than the mean radial flux by a factor ~ a l /3a t ~ 1/6. 
This is to some extent compensated by the fact that a spherical 
core of radius r presents a larger effective cross-section to the 
diffuse flux {4-nr 2 ) than to the radial flux (tit 2 ), i.e. the net dose 
of Lyman continuum radiation received by a spherical core di- 
vides roughly in the ratio 60% radial and 40% diffuse. In order 
to proceed with a spherically symmetric analysis of this phe- 
nomenon, we shall assume that the flux of Lyman continuum 
radiation incident on a core is constant, isotropic, and equal to 
2 (n^^) ■ However, we should be mindful that this is a crude 
representation of the Lyman continuum radiation field incident 
on a real core in a real HII region. We will check retrospec- 
tively that the cores we invoke are much smaller than the HII 
region in which they are imersed (see Section 8). 

4. The photo-erosion rate 

We denote the instantaneous radius of the neutral core by r w , 
and the density of hydrogen nuclei in all forms in the neutral 
core by «i(r), r < r 1F , where r is distance measured from the 
centre of the core. We assume that the ionization front eats into 
the core at speed [- r IF ] < fli , where a x is the isothermal sound- 
speed in the neutral gas. Consequently the ionization front is 
not preceded by a strong shock front, just a weak compression 
wave. This is a reasonable assumption, except near the begin- 
ning, when the ionization front advances very rapidly into the 
core. We can only allow for this by performing numerical sim- 
ulations. However, this initial phase of rapid erosion does not 
last long, and will not greatly affect the overall phenomenol- 
ogy of core erosion or the final protostellar mass. It will simply 
shorten the time interval, f,, before the protostar forms. Thus 
the expressions for the final protostellar mass which we derive 
later (Eans.l40land l4U should be taken as indicative rather than 
precise. 

Conservation of mass across the ionization front requires 



ni(r lF ) [-/if] = n n (r 1F ) [v„(r IF ) - r lF ] , 



(6) 



where «n(r IF ) and Vn(r IF ) are - respectively - the density of hy- 
drogen nuclei in all forms and the radial speed of the ionized 
gas, just outside the neutral core. 

To proceed, we assume that the ionized gas flows away 
from the ionization front at constant speed 



where a u is the isothermal sound-speed in the ionized gas. (In 
reality (Kahn 1969) the flow of ionized gas downstream from 
the ionization front accelerates away from the core due to the 
pressure gradient, but this only changes the result we derive 
below by a numerical factor of order unity.) Hence the density 
in the ionized gas falls off approximately as r~ 2 , 



nn(r) 



«ii(»if) — 
I'd, 



(8) 



until it approaches the background density n at radius r x » 
r lF . Also, substituting from Eqn.Qinto Eqn.|6]gives 

ni(r w ) [-r 1F ] - n n (r w ) «ii • (9) 

Almost all the incident flux of Lyman continuum photons 
is used up balancing recombinations in the outflowing ionized 
gas, so we can put 

;.kv,.„.\ I _2,_w_ „ a*nl(r w )r w 

"Lyc 



2(AC SC ) - f 



a„ n u (r) dr 



(10) 



nnfap) - 



a, r lF 



1/2 



rw - - 



nmirw)an 
rii(rw) 



1 6 Kr)i 



a» r w 



1/2 



All 

niirw) 



(11) 



(12) 



In obtaining the final expression in Eqn. (II 01 . we have assumed 
that r x » r IF . Again this is a reasonable assumption except near 
the beginning. 

5. The first phase (0 < t < u) 

5.1. The initial core 

For the purpose of this exploratory calculation, we assume that 
the pre-existing prestellar core is a singular isothermal sphere. 
Generalizing the treatment to a non-singular (i.e. Bonnor- 
Ebert) isothermal sphere should not change the result greatly, 
but would necessitate a detailed numerical formulation of the 
problem; it would tend to reduce the mass of the final protostar, 
and hence assist in producing low-mass final objects. 

In a singular isothermal sphere, the density of hydrogen nu- 
clei in all forms, n^r) and the mass interior to radius r, M(r), 
are given by 



«i(r) = 



M(r) = 



2nGmr 2 ' 
2 a? r 



(13) 



(14) 



where m = m p /X = 2.4 x 10~ 24 g is the mass associated with 
one hydrogen nucleus when account is taken of other elements 
(in particular helium), m p is the mass of a proton, and X is 
the fractional abundance of hydrogen by mass. If the bound- 
ary density is n Q , the outer radius and total mass are initially 



r = 



(2nGn m) 1/2 ' 

1/2 



[v n (r) - /if] = [vii(rip) - /if] 



(7) 



7rG 3 n Q m 



(15) 



(16) 
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5.2. The inward propagating compression wave 

When this core is first overrun by the HII region, at time t = 0, 
the resulting increase in the external pressure drives a com- 
pression wave into the core, thereby triggering its collapse. We 
assume that the compression wave travels ahead of the ioniza- 
tion front. It leaves in its wake a subsonic - and approximately 
uniform - inward velocity field. The compression wave reaches 
the centre at time 

1 

(17) 



u = r ± = 

fli (2nGn m) 111 



i.e. after one sound-crossing time. We can now justify the as- 
sumption that the time-scale on which the HII region overruns 
the core is much shorter than a sound-crossing time. The time- 
scale for the establishment of the Initial Str0mgren Sphere is of 
order the recombination time-scale, viz. 

^recombination (^*^o) > 



so the condition becomes t KC0 , 



<sc t\ , or 



n » 



2nGm 



3x 1(T 5 cm -3 . 



This is easily satisfied since typically n ~ 1 cm 



6. The second phase (u <t < t 2 ) 

6.1. Formation of the central protostar 

At t = ?i a central protostar forms and subsequently grows in 
mass at a constant rate 



(18) 



We note that this result is true for the collapse of a singular 
isothermal sphere irrespective of whether its collapse is trig- 
gered - as here - from outside-in by the inward propagation of 
a mild compression wave (see Whitworth & Summers, 1985), 
or - as in the Shu (1977) model - from inside-out by a central 
perturbation. 

6.2. The expansion wave 

At the same time that the compression wave converges on the 
centre, an expansion wave is launched outwards at speed oi, rel- 
ative to the gas. If we neglect the small inward velocity already 
acquired by the gas due to the inward propagating compression 
wave, the radius of the expansion wave is given by 



r EW - fl I ( f - f l) - «I f - r ■ 

We now introduce dimensionless variables 
r 

? = -, 

r 

t 

T — — , 

fl 



(19) 

(20) 
(21) 



in terms of which the radius of the expansion wave is given by 



£ew = T - 1 , 



T > 1. 



(22) 



6.3. The core interior to the expansion wave 

Interior to the expansion wave, material flows inwards, to 
feed the constant accretion rate onto the central protostar, and 
quickly approaches a freefall-like profile. Therefore the density 
and velocity can be approximated by 



n[(r) = (3/8nGrn)a'r^r 



2 _-l/2 --3/2 



v'{r) 



(2/3) a^r 



1/2 r -l/2 



t>t 1 , r<r EW . (23) 



There is a degree of freedom in specifying the coefficients of 
nf(r) and vj(r). We have chosen them so that (i) they are consis- 
tent with the constant accretion rate (see Ean.ll8>. and (ii) the 
mass of the central protostar and the mass in the accretion flow 
interior to the expansion wave are equal. Formally the velocity 
changes discontinuously across the expansion wave, from zero 
outside to -2aj/3 inside; in reality the change will be much 
smaller since the inward propagating compression wave during 
the first phase will have set up an inward velocity of compara- 
ble magnitude (e.g. Whitworth & Summers 1985). 

6.4. The ionization front exterior to the expansion wave 
Substituting for «i(r IF ) in Eqn. dl 2l i from Eqn. dl 31 . we obtain 



dr. 



IF 

dt 



1/2 



2nG m an rz 



(24) 



If we now introduce a third dimensionless variable 

1/2 



P = 



° V v Lyc / 'o nGma n 
a, J a] 

(n3«a 6 G 3 m 3 Nt 



v Lyc 



2 7 af n 



an 
2ai 



Eqn. J24l> reduces to 
with solution 

& = [i + PtY 2 ■ 



(25) 



(26) 



(27) 



To fix the constant of integration, we have used the initial con- 
dition, viz.'. at t = 0, r lF = £ IF r = r ; and hence at t = 0, 
£f = 1. 

6.5. The ionization front meets the expansion wave 

The ionization front meets the expansion wave at time U when 
= Jew, or in terms of the dimensionless variables, at r 2 when 
£if = £ew = Equating Eqns. (I27> and i22\ . we obtain 

(i + fjT 2 r 2 = t 2 - 1 . (28) 

With typical values A/L yt = 10 50 s _1 and n = 10 4 cm~ 3 , we 
have (N^) ~ 3 x 1°" cm" 2 s _1 . Putting a x = 0.3 km s^ 1 and 
a u = 10 km s~ 1 , we obtain f3 25 . Hence, if we put 



ft = t 2 — 1 , 



(29) 
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it follows that £ 2 ^ 1, and Eqn.J28> becomes 

[(i+js) +m 2 t2 = i, oo) 

with approximate solution 

£ ~ (l + y?)- 2 ~ 0.002. (31) 

This is the dimensionless radius of the inward propagating ion- 
ization front when it meets the outward propagating expansion 
wave. 

7. The third phase (t 2 <t < t 3 ) 

7.1. The ionization front interior to the expansion wave 

Once the inward propagating ionization front passes the out- 
ward propagating expansion wave, the ionization front encoun- 
ters infalling gas with density and velocity given by Eqn. l23l 
Eqn. [9] therefore has to be adjusted to account for the inward 
motion of the neutral gas: 

n[{r w ) [-/"if + v[(r w )} «n(r IF ) a u . (32) 

Substituting for n[(j\p) and v,'(ri F ) from Ean.l23land for «n(rm) 
from Eqn. II II we obtain 

= - 2 Jl^l[^ + 4^]. (33) 



The solution is 



(T-l) 



3/2 



■3/2 



80 



I'n 



1 + 4/3£ 



3/2 



1 + 



(34) 



where we have fixed the constant of integration by requiring 
that r IF = r EW = r 2 at t 2 \ or equivalently £ IF = £ FW = £ 2 = t 2 - 1 
at t 2 . 

To follow the ionization front into the third phase, we pick 
ever smaller values of £ IF (< ij z ) and calculate r from Eqn. 1341 
The radial velocity of the ionization front is then given by r 1F = 
fli^F, where £ IF can be calculated from Eqn. [33] The specific 
kinetic energy of the newly ionized gas, relative to the central 
protostar is therefore 



£fCIN 



1 2 



?IF H 



(35) 



Taking into account the mass of the protostar, M* = 
1), plus the mass in the accretion flow interior to 
the ionization front, M, = G~ l c%r (j - ly 1 ' 2 * 3 ! 2 
binding energy of this material is 



G l a 2 r Q (j- 



^ft^, the specific 



Cbind 



G{M„ + Mi) 



(T-l) 

ft* 



,1/2 



(r- 1)1/2 



(36) 



We shall assume that ionization ceases to erode the core once 
the specific kinetic energy, e K iN (Eqn.l35>. is less than the spe- 
cific binding energy, e BIND (Eqn.l36>. The mass M 3 = M, + M, 
interior to the ionization front at this stage is the final mass of 
the protostar. This occurs at time t 3 = t 3 ?i, where 



6f(t 3 ) + 



(r 3 - 1) 



ft 



1/2 



ft (t 3 - \yi 2 ' 

and the final mass of the protostar is then given by 



(37) 



(38) 



8. Results 

The procedure for obtaining the final protostellar mass M 3 is 
therefore as follows, (i) Chose (N Lyc ,ai,n ) and calculate [} 
from Eqn|25] (ii) Solve Eqn. H5|for ft. (iii) Solve Eqn. E] 
for^ IF (r), (t > 1 + ^ 2 )- (iv) Identify t 3 satisfying Eqn. 1371 (v) 
Compute M 3 according to Eqn. 1381 

Figure 1 shows, for three representative values of // Lyc [i.e. 
(a) 10 50 s _1 , (b) 10 51 s -1 , and (c) 10 5 V] loci of contantM 3 = 
O.OO1M , O.O12M , and O.O78M on the (ai,fogio[«o])- 
plane. These values of M 3 are chosen because they represent, 
respectively, one Jupiter mass, the deuterium-burning limit, 
and the hydrogen-burning limit. Notionally objects in the mass 
range O.O12M $ M 3 $ O.O78M are referred to as brown 
dwarves, whilst those in the mass range M 3 ^ O.O12M are 
referred to as planetary-mass objects. We see that for a broad 
range of values of A/Lyc, «o and a,\, there exists the possibility to 
form free-floating brown dwarves and planetary-mass objects, 
provided suitable pre-existing cores are available to be photo- 
eroded. 

In estimating the radiation field incident on a core, we have 
assumed that the volume occupied by the core is much smaller 
than the volume of the HII region around it (see Section 3). It is 
therefore necessary to check that this is the case. If we require 
the ratio between the initial volume of a core and the volume 
of the whole HII region to be less than some fraction 77, then 
substituting from Eqns.lT1andll5lwe obtain the condition 



«o < «„ 



n 

M = 5 



3M 



Lyc?7 



a* 



G m 

„2 



(39) 



« max (77) is plotted on Figure 1 for 77 = 0.001 (dashed curve) 
and rj = 0.01 (dash-dot curve). Evidently this is a significant 
constraint for VV Ly( . ~ 10 50 s" 1 , but for higher values of A/Lyc it 
is less important. 

Since the protostar does not grow hugely during the Third 
Phase, we can approximate the results by adopting the mass at 
the end of the Second Phase. If in addition we assume /3 » 1, 
so that from Eqn. [2] ft - P~ 2 , we can write 



M 3 fM t t^ 2 
f2 3 a, af 



1/3 



0.0 \M ( 



'\0.3kms-V Uo 5 ^- 1 / UOWV { ' 



Here / ^ 1.7 is a factor to account for the additional mat- 
ter which accretes onto the protostar after the ionization fron 
meets the expansion wave. We can understand this factor as fol- 
lows. In a collapsing singular isothermal sphere, the infalling 
mass interior to the expansion wave is equal to the mass al- 
ready in the central protostar. Thus, in the problem considered 
here, when the ionization front meets the expansion wave at 
t = t 3 , the infalling mass interior to the ionization front equals 
the mass already in the protostar. It follows that if no further 
mass acretes onto the protostar, / = 1; and if all the infalling 
mass interior to the ionization front at t = t 3 accretes onto the 
protostar, then / = 2. The empirical value / = 1.7 suggests 




Fig. 1. For (a) N Lyc = 10 50 s -1 , (b) N Lyc = 10 51 s -1 , and (c) N Lyc = 10 52 s -1 , the solid lines show the loci on the (a u {og m [n ])- 
plane where free-floating objects of mass 0.001 M Q (one Jupiter mass), 0.012 M (the deuterium-burning limit) and 0.078 M Q (the 
hydrogen-burning limit) form due to photo-erosion of pre-existing cores. The dotted lines show the predictions of the approximate 
Ean.1401 The dashed (dash-dot) line is the locus below which the volume of the initial core is less than 0.001 (0.01) of the volume 
of the whole HII region. 



that most of the infalling mass interior to the ionization front at 
t — f 3 is accreted onto the protostar. 

Ean.l40lis plotted as a dotted line on Figure 1. We see that 
it is a very good approximation to the numerical results. Given 
the approximate nature of the whole analysis, Eqn.@U|can be 
viewed as our best estimate of the final protostellar mass. 

9. Discussion 

This exploratory analysis suggests that it is possible for a rel- 
atively massive prestellar core which suddenly finds itself im- 
mersed in an HII region to be stripped down to the mass of a 
brown dwarf, or even a planet, by photo-erosion, before it can 
collapse to form a star. To improve on this analysis, it is proba- 
bly necessary to resort to numerical modelling. One could then, 
in principle, (i) use a more realistic configuration for the pre- 
existing prestellar core; (ii) include the development of the HII 



region as it overruns the core; (iii) treat properly the early phase 
when the ionization front advances very rapidly into the outer 
layers of the core, and the compression wave which is driven 
into the core; (iv) attempt to capture effects due to departures 
from spherical symmetry, in particular the radiative transfer 
aspects; (v) include self-gravity explicitly in the equations of 
motion. We suspect that the photo-erosion mechanism may be 
somewhat more effective than our analysis suggests, because 
firstly, with a prestellar core which is not as centrally condensed 
as a singular isothermal sphere (e.g. Ward-Thompson, Motte & 
Andre, 1999; Bacmann et al. 2000), the protostar will not form 
so quickly, giving the ionization front more time to erode; and 
secondly, the ionization front will penetrate at higher speed into 
the less dense central regions. 
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We can express the final protostellar mass M 3 in terms of 
the initial core mass M , by combining Eans. I40lll6l andmto 
obtain 



2fa,GM 2 
3 a x a\ R m i 



There are two important deductions which we can make from 
this equation. 

First, the ratio of the final protostellar mass to the initial 
core mass is given by 



Ml „ lQ -A( M o\( R wa 



M 



Mr. 



pc 



(42) 



The photo-erosion process is therefore very effective, in the 
sense that it strips off most of the mass of the initial core. 
This makes it a rather wasteful way of creating low-mass pro- 
tostars, in the sense that it requires a relatively massive initial 
core to create a single final low-mass protostellar object. In ad- 
dition, we infer that any intermediate-mass protostars which 
have formed in the vicinity of a group of OB stars must already 
have been well on the way to formation before the OB stars 
switched on their ionizing radiation; otherwise these protostars 
would have been stripped down to extremely low mass. 

Second, if the number of cores in the mass interval 
(M Q , M Q + dM ) (i.e. the core mass function) has the form 



N Mo dM oc M Q y dM , 



(43) 



then the initial mass function for brown dwarves and planetary- 
mass objects formed by photo-erosion has the form 



N M3 dM 3 oc M 3 (r+1)/2 dM 3 . 



(44) 



Motte et al. (1998) find y ^ 1.5 for M < O.5M and y 
2.5 for M > O.5M for the cores in p Ophiuchus, and similar 
results are reported for the cores in Ophiuchus by Johnstone et 
al. (2000), for the cores in Serpens by Testi & Sargent (1998), 
and for the cores in the Orion B molecular cloud by Motte et 
al. (2001) and Johnstone et al. (2001). Hence the exponent in 
the initial mass function for brown dwarves and planetary-mass 
objetcs formed by photo-erosion should lie in the range - 1 .25 
to - 1 .75 , with the smaller value (-1 .75) obtaining for higher- 
mass cores. As the HII region expands and engulfs additional 
cores, R HU increases, and therefore for fixed M 3 the progenitor 
cores have larger masses M . Consequently, as one moves away 
from the stars exciting the HII region, the density of brown 
dwarves and planetary-mass objects formed by photo-erosion 
should decrease, and their initial mass function should become 
slightly steeper. 

We have omitted from our analysis the possible effect of 
a stellar wind impinging on a core (A Burkert, private com- 
munication). A simple estimate suggests that a stellar wind 
would need to be very powerful to have a significant effect. 
Suppose that the central OB star blows a wind with mass-loss 
rate M w i„d ~ 10~ 4 M Q year" 1 , and speed v w ; n d ~ lOOOkms" 1 . 
The resultant ram pressure acting on a core at radius R is 



wind — 



^^wind ^wind 

4ttR 2 



(45) 



We can compare this with the ram pressure of the ionized gas 
flowing off the surface of a core which is exposed to the direct 
flux of ionizing radiation: 



(41) P 



IF 



7V L y C m a n 



AnR 1 ' 

Substituting typical values, we obtain 



(46) 



^^wind ^wind 



■flF A^Lyc m an 



= 6xlQ- 



lO- 4 M year 



/ I'wind \ 

1 1000 km s- 1 / 



Nhyc 



[5 x lO^s" 1 J 



(47) 



This implies that the wind can have little effect. In reality, most 
of the ionizing flux from the central star is spent maintain- 
ing ionization against recombination in the flow of ionized gas 
from the core, rather than producing new ionizations. Therefore 
Eqn. J46l > significantly overestimates the ram pressure of the 
ionized gas flowing off the core, but probably not by two orders 
of magnitude. The effect of a stellar wind will be to modify the 
rate at which the ionized gas flowing off the core can disperse 
(e.g. Richling & Yorke, 1998; Hollenbach, Yorke & Johnstone, 
2000). A bow shock will form on the side of the core facing 
into the wind (the windward side), and the ionized gas flow- 
ing off this side of the core will be swept back in the direction 
of the wind. This will have the effect of increasing the density, 
and hence the recombination rate, in the ionized gas flowing off 
the core, thereby reducing the rate of erosion somewhat. On the 
leeward side of the core the effect of the wind will be minimal. 

10. Conclusions 

The analysis we have presented in this paper is based on several 
assumptions and approximations, which, whilst reasonable, are 
idealizations of what is likely to occur in nature. With this pro- 
viso, we infer that, in the central regions of large clusters, free- 
floating brown dwarves and planetary-mass objects may form 
when more massive pre-existing cores are overrun by the HII 
regions excited by newly-formed massive stars. The sudden 
ionization of the surroundings of a core drives a compression 
wave into the core, creating at its centre a protostar which then 
grows by accretion. At the same time an ionization front starts 
eating into the core, thereby removing its outer layers. The fi- 
nal mass of the protostar is determined by a competition be- 
tween the rate at which it can accrete the infalling envelope, 
and the rate at which the ionization front can erode the enve- 
lope. The simple analysis presented here suggests that in large, 
dense clusters this will result in the formation of free-floating 
brown dwarves and/or planetary-mass objects. The process is 
both robust, in the sense that it operates over a wide range of 
conditions, and also very effective, in the sense that it strips 
off most of the mass of the initial core. This means that it is a 
wasteful way of creating low-mass protostars, because a rela- 
tively massive initial core is required to produce a single final 
brown dwarf or planetary-mass object. It also means that any 
intermediate-mass protostars which have formed in the vicin- 
ity of a group of OB stars must already have been well on the 
way to formation before the OB stars switched on their ionizing 
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radiation; otherwise these protostars would have been stripped 
down to extremely low mass. 
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